The Role of Electron-electron Interactions in Graphene ARPES Spectra 
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We report on a theoretical study of the influence of electron-electron interactions on ARPES 
spectra in graphene that is based on the random-phase-approximation and on graphene's massless 
Dirac equation continuum model. We find that level repulsion between quasiparticle and plasmaron 
resonances gives rise to a gap-like feature at small k. ARPES spectra are sensitive to the electron- 
electron interaction coupling strength Qgr and might enable an experimental determination of this 
material parameter. 
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Introduction — In quasi two-dimensional (2D) crystals 
ARPES is a powerful probe of band structure and of 
Coulomb and phonon-mediated electron-electron interac- 
tions. A number of recent experiments [l[ have reported 
ARPES spectra for single-layer and bilayer graphene sys- 
tems grown epitaxially Q on the surface of SiC. These 
recently isolated 3] 2D crystals have attracted consider- 
able attention because of unusual properties 0, Q that 
follow from chiral band states, notably unusual quantum 
Hall effects and because of potential for appHcations. 
Experiment appears to establish [l[ that both electron- 
phonon and electron-electron dressings have an influence 
on graphene's ARPES spectra. Recent theoretical stud- 
ies 0, H, @] provide a solid basis for interpreting the 
electron-phonon contributions; in this Letter we report 
on a theoretical study of the influence of electron-electron 
interactions. Our analysis is based on the massless Dirac 
model which describes the tt electron states within 
~ 2eV of the Fermi energy of a neutral graphene sheet 
and on the random phase approximation (RPA) for the 
self-energy. 

The self-energy in a system of fermions is naturally 
separated into an exchange contribution due to interac- 
tions with occupied states in the static Fermi sea, and a 
correlation contribution due to quantum fluctuations of 
the Fermi sea [ll| . Graphene differs from the widely 
studied 2D systems in semiconductor quantum wells be- 
cause its quasiparticles are chiral and because it is gap- 
less and therefore has interband quantum fluctuations on 
the Fermi energy scale. In graphene band eigenstate chi- 
rality endows exchange interactions with a new source of 
momentum dependence which renormalizes the quasipar- 
ticle velocity and strongly influences the compressibility 
and the spin-susceptibility l^, 3, 15 1. 

Our results [l^ for the ARPES spectra of a n-doped 
graphene sheet with agr = 2 are summarized in Fig. [1] 
Here ctgr — ge^/ehv is the interaction coupling strength, 
9 = 5spin<?vaiicy = 4 is a band degeneracy factor, v is the 
Fermi velocity, and e depends on the dielectric environ- 
ment of the graphene layer (agr ~ 2 is a typical value 
thought to apply to graphene sheets on the surface of a 
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FIG. 1: (Color online) Spectral function A{k,u!) of a n- 
doped graphene sheet for with wavevector k in units of Fermi 
wavevector fcp and energy u in units of and measured from 
the Fermi energy hvkp. The ARPES spectra for each k is 
the portion of the spectral function with u < 0. The k- 
dependence is represented in this figure by results for ten 
discrete k G [0.0,0.9] values separated by 0.1. It follows from 
particle-hole symmetry that the spectral function of a p-doped 
graphene sheet is specified by the tj > spectral weight. 



Si02 substrate). The most notable aspect of these re- 
sults is the appearance near fc = of strong plasmaron 
peaks in addition to quasiparticle peaks, which give rise 
to a gap-like structure in the overall spectrum. In the 
following paragraphs we explain the physics behind this 
figure. 

Doped Dirac Sea Charge Fluctuations — The massless 
Dirac band Hamiltonian of graphene is 



H = VT ((Tl pi -I- 0-2 P2) , 



(1) 



where r = ±1 for the inequivalent K and K' valleys 
at which tt and tt* bands touch, pi is an envelope func- 
tion momentum operator, and ai is a Pauli matrix which 
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FIG. 2: (Color online) Left panel: —Qm[e~^{q,u!)] as a func- 
tion of q/kp and uj/ef for Qgr — 2. The red solid line is the 
plasmon dispersion relation, u)pi{q ^ 0) = £f ^/a^rq/J^kp) . 
Right panel: ,uj)] as a function of q/kp and 

uj/ef- The left and right panels become identical in the non- 
interacting agr — > limit. 



acts on the sublattice pseudospin degree-of-freedom. The 
low-energy valence band states have pseudospin aUgned 
with momentum, while the high energy conduction band 
states, split by 2u|p|, are anti-aligned. In Fig. [2] we 
compare the particle-hole excitation IJ, ll7| spectra of 
non-interacting and interacting 2D doped Dirac systems. 
The non-interacting particle-hole continuum is repre- 
sented here by the imaginary part of graphene's Lind- 
hard function, 5rn[x(°)(<7, w)], which weighs transitions 
by the strength of the density fluctuation to which they 
give rise; transitions between states with opposite pseu- 
dospin orientation therefore have zero weight. More 
generally the band-chirality related density-fluctuation 
weighting factor (called the chirality factor below), which 
plays a key role in the physics of the ARPES spectra, is 
[1 ± cos(0fc.fc+q)]/2 with the plus sign applying for intra- 
band transitions and the minus sign applying for inter- 



band transitions, and 9k,k+q equal to the angle between 
the initial state (fc) and final state {k + q) momenta. 
The weight is therefore high for intraband (interband) 
transitions when k and k + q are in the same (opposite) 
direction. The main features in Fig. [2] that are impor- 
tant for ARPES spectra are the l/\/vq — divergence 
which occurs near the upper limit of the q < kp in- 
traband particle-hole continuum and the relatively weak 
weight at the lower limit of the q < kp inter-band 
particle-hole continuum. The divergence at the intra- 
band particle-hole spectrum contrasts with the singular 
but finite V^max — ^ behavior at the upper end of the 
particle- hole continuum in an ordinary electron gas. The 
difference follows from the linear quasiparticle dispersion 
which places the maximum intraband particle-hole exci- 
tation energy at vq for all k in the Fermi sea. 

In the RPA, quasiparticles interact with Coulomb- 
coupled particle-hole excitations. Because the bare 
particle-hole excitations are more sharply bunched in en- 
ergy. Coulomb coupling leads to plasmon excitations that 
are sharply defined out to larger wavevectors than in the 
ordinary electron gas and steal more spectral weight from 
the particle-hole continuum. As seen in Fig. [21 the plas- 
mon excitation of the Dirac sea remains remarkably well 
defined even when it enters the interband particle-hole 
continuum. The persistence occurs because transitions 
near the bottom of the interband particle-hole contin- 
uum have nearly parallel k and fe -f- g and therefore little 
charge-fluctuation weight. Interactions between quasi- 
particles and plasmons are stronger in the 2D massless 
Dirac system than in an ordinary non-relativistic 2D sys- 
tem. 

Dirac Quasiparticle Decay — In Fig. [3] we plot our results 
for the imaginary part of the Dirac sea self-energy for a 
series of fe values (these results are for a n-doped agr = 2 
system). In the RPA 



j^-^ Vq 3TO[e-^(q,w-^^/(fc + q)) 



1 + SS' cos {0k.k+q) 



[e(^-es'(fc + g)) 



e(-6'(fc + q))] 

(2) 



where s,s' = ±1 are band indices, Vq = 27re^/(eg) is the 
2D Coulomb interaction, e{q,uj) — 1 — Vqx'^^Hqjto) is the 
RPA dielectric function, and 8(x) is the Heaviside step 
function. The two factors in square brackets on the right- 
hand-side of Eq. ([2]) express respectively the influence 
of chirality and Fermi statistics on the decay process. 
Note that depends on the band-index s only through 
the chirality factor. For a; > and fixed q, the RPA 
decay process represents scattering of an electron from 
momentum fc and energy w to fe -f q and ^s'{k + q), 
with all energies in Eq. ([2]) measured from the Fermi 



energy. Since the Pauli exclusion principle requires that 
the final state is unoccupied, it must lie in the conduction 
band, i.e. s' = +1. Furthermore since the Fermi sea is 
initially in its ground state, the quasiparticle must lower 
its energy, i.e. £,s' < ^ - electrons decay by going down in 
energy. For a; < 0, the self-energy expresses the decay of 
holes inside the Fermi sea, which scatter to a final state, 
by exciting the Fermi sea. In this case the final state 
must be occupied so both band indices are allowed for 
s' , and energy conservation requires that holes decay by 
moving up in energy. Since photoemission measures the 
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FIG. 3: (Color online) The absolute value |Sfm[Es(fe,aj)]| of 
the imaginary part of the RPA quasiparticle self-energy (in 
units of ep) of an n-doped system as a function of energy uj 
for fc = 0, 0.25, 0.50, and 0.75 and for Qgr = 2. 



properties of holes produced in the Fermi sea by photo 
ejection, only < is relevant for this experimental 
probe. Note however that because of the particle-hole 
symmetry properties of the Dirac model I]s(fc,cj) in an 
n-doped system is identical to I]_s(fe,— tj) in a system 
with the opposite doping. Our results for 3to[Ss(A;, w)] 
for w > 0, therefore specify the ARPES spectra of p- 
doped graphene. 

We remark that because interaction and band ener- 
gies in graphene's Dirac model both scale inversely with 
length, 5m[Es(fe, w)] = vk-pF{u)/vk-p,k/k-p). For large 
|a;|, F{x,y) -Tra^ i(agr)|x|/(64g), where 



("gr) 



dx 



1 



a;2/256 



[^(0) = 4/3, l{2) ~ 0.655124]. This implies that for \uj\ > 
wfcp, the decay rate in a doped system (5TO[Ss(fc, oj)]) ap- 
proaches that of an undoped system (the Fermi energy 
£f — vk-p is used as the energy unit and fcp as the unit 
of wavevector in all plots and in the remaining sections 
of this Letter). As we will see, however, the doped sys- 
tem properties are quite different from those of an un- 
doped system up to energy scales several times larger 
than the Fermi energy, particularly so near the Dirac 
[k = 0) point. 

In explaining the spectra plotted in Fig. [3] we start 
with the fc = case (9m[E+(0,a;)] = 3to[E_(0, t^)]) for 
which the final state energy is (q) — s'q—1. For lu > 
the final state must be unoccupied so that s' = -1-1; g is 
restricted to those values larger than 1 for which the the 
Dirac sea excitation energy il{q) = lu + I — q is positive. 
Comparing with Fig.[2]we see that 5m[I]_|_(0, lu)] vanishes 
like Lu'^ for w — s- 0, a universal property of normal Fermi 



liquids. In the RPA 



3to[S^(0, lu 



0)] = - 



x/3 



gr 



8g (l + agr)2 



(4) 



The sharp increase in !3to[S+(0, cj)] which occurs at 
LU 1.2 reflects the onset of plasmon emission. Note 
that these plasmons at g > 1 remain well defined exci- 
tations well into the interband particle-hole continuum. 
For u! < 0, relevant to ARPES in n-doped graphene, 
both conduction and valence band final states occur and 
transitions are allowed if the transition energy ^{q) = 
\lu\ — 1 + s' q is positive and the final hole state is occu- 
pied, i.e. q < 1 for s' = +. Note that for final states in 
the conduction band, the transition energy is close to the 
plasmon excitation energy over a wide range of q values. 
This property leads to much stronger plasmon features 
for LU < than for lu > 0. 5m[I]+(0,a;)] diverges for 
LU ~ —1.3, the value at which fl{q) for s' = +1 is tangent 
to the plasmon dispersion line illustrated in Fig. [21 In 
this case, the hole scatters into a resonance consisting of 
a quasiparticle strongly coupled to plasmon excitations, 
a plasmaron [3, [l^. At finite fc, the conduction and va- 
lence band Qto[Es(A;, w)] plasmaron peaks broaden and 
separate, because of the dependence on scattering angle 
of ^s' (k+q) and because of chirality factors which empha- 
size k and q in nearly parallel directions for conduction 
band states and k and q in nearly opposite directions for 
valence band states. Because of the chirality factor the 
conduction band plasmaron moves up in energy approx- 
imately as vk while the valence band plasmaron moves 
down. 

ARPES Spectra — We are now prepared to discuss the 
ARPES spectra in Fig. [TJ ARPES measures the wavevec- 
tor dependent electron spectral function 



(3) Asik,Lu) 



|3ml]s(fe, Lu)\ 



1 

^[iu-^sik)- 3fieEs(fc,cj)]2 + ['^m^s{k,Lu) 



(5) 

Near the Fermi energy the spectral function consists of a 
narrow Lorentzian centered at the energy E which solves 
the Dyson equation for the quasiparticle energy: 



E ^Uk) + 'fie^sik,E). 



(6) 



For small fc we see in Fig. |4] that there are two solutions 
to this equation, which are shifted from the bare quasi- 
particle energy and the plasmaron energy because of the 
strong energy-dependence of SReEs near the StoEs's plas- 
maron peak. These shifts can be understood as following 
from level repulsion between a bare particle resonance 
and a plasmaron resonance. We also note in Fig. [5] that 
SReEs has a negative contribution which is present at the 
Fermi energy and persists over a wide regime of energy. 
This contribution is due to exchange and correlation in- 
teractions of quasiparticles near the Fermi energy with 



the negative energy sea. As explained 



pre- 
viously, this effect produces a nearly rigid shift in the 
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FIG. 4; (Color online) »e[E+ (fc, u)], 9m[E+ (fe, a;)], and spec- 
tral function A+{k,Lj) for k — 0.25 and fc = 0.75. The band 
energy and 5ReE+ are measured from the band ep and interac- 
tion [E+(fcF,a; = 0)] contributions to the chemical potential. 



band energies which is increasingly negative further be- 
low the Fermi energy, increasing the band dispersion and 
the quasiparticle velocity. In Fig. 3] we show that for 
larger fc, the plasmaron feature in 3mEs is broadened 
sufficiently to remove the plasmaron solution of the quasi- 
particle Dyson equation. In this limit, the spectral func- 
tion has only a weak plasmon satellite on top of the main 
quasiparticle peak. 

The full ARPES spectral in Fig. [1] can be understood 
by applying considerations like those explained above to 
both bands and over the full range of k. Because the 
negative exchange-correlation shift in the quasiparticle 
energy gets larger as A: —> in the conduction band, 
the energy splitting between the sharpening plasmaron 
resonance and the band-like quasiparticle decreases and 
their coupling increases. In addition, because the slope 
difference between E — £,s{k) and dieY^sik, E) is smaller 
at the plasmaron Dyson equation solution than at the 
quasiparticle Dyson equation solution, the plasmaron as- 
sumes the larger fraction of the spectral weight. The 
detailed behavior near fc = is sensitive to the value of 
Qfgr. For agi- = 2, the avoided crossing between quasipar- 
ticle and plasmaron solutions occurs near k = 0. Because 
the higher energy valence band peak has little spectral 
weight when it crosses the Fermi energy, the overall spec- 
trum can appear to have an energy gap. 

Because electron-electron interaction effects occur on 
the vkp energy scale they can be separated from electron- 
phonon interaction effects experimentally by varying car- 
rier density. Although the RPA is not exact, it pro- 
vides a good starting point for interpreting the influence 
of electron-electron interactions on the ARPES spectra 
of graphene sheets. As we have shown, the downward 
shift of valence band states compared to conduction band 
states and plasmaron-quasiparticle level repulsion and 
weight transfer effects must be accounted for to interpret 
ARPES data. Features in the data which might other- 
wise be interpreted as extrinsic effects related to disorder 
or interactions with substrates can have interesting and 
subtle intrinsic origins. Close comparison between the- 
ory and experiment and further materials progress will 



augment the power of ARPES experiments. 
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